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Internal waves are believed to be of primary importance as they affect ocean mixing 
and energy transport. Several processes can lead to the breaking of internal waves and 
they usually involve non linear interactions between waves. In this work, we study ex- 
perimentally the parametric subharmonic instability (PSI), which provides an efficient 
mechanism to transfer energy from large to smaller scales. It corresponds to the desta- 
bilization of a primary plane wave and the spontaneous emission of two secondary 
waves, of lower frequencies and different wave vectors. Using a time-frequency analy- 
sis, we observe the time evolution of the secondary waves, thus measuring the growth 
rate of the instability. In addition, a Hilbert transform method allows the measurement 
of the different wave vectors. We compare these measurements with theoretical predic- 
tions, and study the dependence of the instability with primary wave frequency and 
amplitude, revealing a possible effect of the confinement due to the finite size of the 
beam, on the selection of the unstable mode. 



1. Introduction 

An essential ingredient of the thermohaline circulation is the mechanism by which the 
denser water that was produced in the high latitude regions (colder and saltier), once 
it has flowed down the continental slopes towards the abyssal ocean, can come back 
to the surface to close the loop. This process involves an energy input to provide the 
gain of potential energy necessary to lift this denser water. It is believed that turbulent 



mixing, generated by wind and tides, is the mechanism that performs this task ( Munk 



1966 Munk & Wunsch||1998). One possible way is via the breaking of internal gravity 



waves |Staquet & Sommeria 2002| , ubiquitous in the ocean, allowing a transfer of energy 



from large scales to small scales, where this energy is partly dissipated in heat and partly 
converted in potential energy through diapycnal mixing. 

The detailed mechanisms for energy dissipation of internal gravity waves is still de- 
bated. Several mechanisms have been proposed: 

• parametric subharmonic instability (|MacKinnon & Winters 2005 Alford et al. 2007| , 



reflection on sloping boundaries (Dauxois & Young 1999 |, ~ 
scattering by mesoscale structures ( Rainville & Pink el|2006b|, 
small scale bathymetry ( Kunze & Smith 2004,: Johnston et al.|[2003 Peacock et al. 
2009) . ' ' 



The importance of these four possible dissipative processes has to be estimated and 
compared precisely. It is likely that a combination of them might be the correct answer 
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but the usual physicists' approach, aiming at separating the different processes one by 
one, is presumably appropriate in a first stage. 

Parametric subharmonic instability (PSI) is the resonant mechanism by which a pri- 
mary wave is unstable to infinitesimal perturbations, transferring energy through the 
quadratic nonlinearity of the Navier-Stokes equation to two secondary waves, satisfy- 
ing temporal and spatial resonant conditions. More precisely, PSI is the class of resonant 
wave-wave interactions wherein energy is transferred from large scale to smaller scales 
and where the frequencies of the secondary waves are near half the primary frequency. 
Studying the particular case of a primary plane wave is very interesting since a plane 
wave is solution of the full rnviscid nonlinear equation for any amplitude. However, as 
will be discussed below, this solution becomes unstable above a given threshold, which 
can be computed analytically (McEwan et al. 1972 1. 



Initially, this instability was considered only for the gravity or capillary waves (IMc- 
( 1965, )). However, Hasselmann ( |1967[ | proved that this instability could be ob- 
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served in many physical systems, and in particular for internal waves. Following this 
initial report, several theoretical studies have been developed for this phenomenon, 
deriving in particular the expression for the growth rate in presence (or absence) of 
viscosity ( Thorpe|[l968 



McEwan 



1971 McEwan et al. 



1972 



McEwan & Plumb||1977 1. 



At the same time, the first experimental observations of the instability have been re- 



ported (.McEwan|1971{ [McEwan et al.|1972} [McEwan & Plumb|1977-,Benielli & Somme- 



ria |1998) . These first observations allowed the determination of the amplitude threshold 
of the instability and the temporal evolution of the secondary waves amplitude. With 
the development of more powerful computers, numerical simulations have been devel- 
oped on this subject (|Bouruet-Aubertot et al. 1995 Camevale et al.||2001 Koudella & 
|Staquet'2006 ) to estimate the energy transfer between the different scales and to deter- 
mine different scaling laws. Finally, only very recently new experiments using a vertical 
mode-1 wave have been performed by Joubaud et al. ( 2012| , allowing precise measure- 
ment of the frequencies and the wave vectors of the secondary waves and to compare 
theses measurements to the viscous theory. Note that a study of PSI in the very similar 
context of rnertial waves in a rotating fluid has recently been reported by Bordes et al. 
( |20T2) . 

In this paper, we present a comprehensive study of two-dimensional parametric sub- 
harmonic instability. The results of experiments are compared with theoretical predic- 
tions. The paper is organized as follows. The experimental observations for plane waves 
are presented in §|2j followed by the presentation of the theoretical derivation of thresh- 
old, resonant conditions and growth rates in §|3] A careful comparison between experi- 
ment and theory is discussed in §|4] After having studied the dependence with primary 
wave frequency and amplitude in § |5| we present our conclusions and draw some per- 
spectives . 



2. Propagation of plane waves: experimental observations 

2.1. Experimental setup 

A tank, 160 cm long and 17 cm wide, is filled with linearly stratified salt water with 
constant buoyancy frequency N using the standard double bucket method. An internal 
wave is generated using a wave generator similar to the one employed in previous ex- 



periments, described in Gostiaux et al. (2007) and characterized in Mercier et al. (2010 >. 
The generator is composed of stacked plates set in motion thanks to eccentered cams ro- 
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Figure 1: (color online) Sketch of the experimental setup showing the wave generator 
lying horizontally at the top of the wave tank. The dashed lines delimit the expected 
domain of propagation of the left to right propagating wave beam, once it has been 
emitted from the generator. For the sake of clarity, the direction of propagation of the 
phase and group velocities are also indicated. The tilted grey rectangle corresponds to 
the analysis area for the time-frequency study of section |4!2| The plot on the right shows 
the experimentally measured modification of the density due to salt, Ap = p — 1000, as 
a function of the water depth z. Points correspond to experimental measurements while 
the straight line is a linear fit. 



fating inside the plates at constant frequency around a shaft. The motion of these plates 
imposes a moving boundary condition generating the desired wave. 

The main difference with the generator described previously is that the present ver- 
sion allows the eccentricity of each cam to be varied. The same set of cams can there- 
fore be used to generate different profiles, as well as identical profiles with different 
amplitudes. In the experiments presented here, the wave generated is a plane wave of 
wavelength A=72 mm. In this case, the generator is placed horizontally as shown in 
Fig. [T| The motion of the plates is thus defined by a vertical velocity w{x,z = 0,t) = 
acoQ cos(a;of ~ 27rx / A), ci^o being the excitation frequency, a the amplitude of oscillation 
of the plates and H the water depth. Note that to avoid spurious emission of internal 
waves on the extremities of the moving region, the amplitude of the plates is constant 
over two wavelengths in the central region, while one half-wavelength with a smooth 
decrease of the amplitude is added on each sides. 

The motion of the fluid is captured by the synthetic schlieren technique using a dotted 



image behind the tank ( Dalziel et al.|2000^ . A camera is used to acquire images of this 
background at 1.875 frames per second. The CIVx algorithm (j Fincham & Delerce|2000) 
computes the cross-correlation between the real-time and the t = background images, 
when the fluid is at rest. This algorithm gives the variation of the horizontal, px {x, z, t) = 
dx{p{x,z,t) — po{z)), and vertical, pz{x,z,t) = dz{p{x,z,t) — po{z)), density gradients, 
where p{x,z, t) and |0o(z) are the instantaneous and initial fluid densities. 

2.2. Direct observation 

In an experimental configuration where the tank is filled with a linear stratification pro- 
ducing a buoyancy frequency N = 0.91 rad/ s, a plane wave is generated at a frequency 
coq/N = 0.74, with an amplitude of the plates motion of a = 0.5 cm. Figure|2|a) shows 
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Figure 2: (color online) Snapshots of the vertical density gradient field for t = IOTq (a) 
and t — 50To (b) where Tq = In I <xiq is the primary wave period. The wave is propa- 
gating from left to right. On the left panel, the direction of the phase velocity ijt and 
the group velocity are indicated. Note that the shade scale (color scale online) is the 
same in both panels, buoyancy frequency is N = 0.91 rad/s, the wave frequency is 
ojqIN — 0.74 and the motion amplitude of the plates of the generator is 0.5 cm. 



a snapshot of the density gradient field, obtained 10 oscillating periods after the wave 
generator was started. One can observe a plane wave extending over two wavelengths, 
propagating from the top left corner to the bottom right corner of the field of view, while 
the phase propagates from bottom left to top right. Note that the wave has reached a 
steady state in the visualization window. Forty oscillating periods later, a strong per- 
turbation of the plane wave can be observed, as shown in Fig. |2|b) (see the movie in 
the online supplementary material). Smaller scale patterns have formed over the whole 
area that was initially occupied by the plane wave beam; interestingly, these patterns 
extend even slightly outside of this area, in the top left corner. 

2.3. Analysis 

The measured density gradient fields are analyzed using a time-frequency representa- 
tion ( |Flandrin|1999 j calculated at each spatial point 
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where r stands for x or z and h is a smoothing Hamming window of energy unity. A 
large (resp. small) window provides good frequency (resp. time) resolution. To increase 
the signal to noise ratio, the data is averaged on the entire area of observation. In the 
following, we will consider only the analysis of the vertical density gradient field, but 
the results are similar for the horizontal one. 

Figure |3][a) shows the time-frequency spectrum for the experiment corresponding to 
the snapshots presented in Fig. |2] One can clearly observe that initially, only the fre- 
quency coq/N = 0.74 is present: it corresponds to the wave produced by the generator, 
which we will call the primary wave. After about 10 oscillation periods, one notices the 
growth of two secondary waves, with the frequencies coi/N = 0.50 and 0^2/ AT = 0.24. In 




Figure 3: (color online) (a) Time frequency spectrum Sz{co,t) of the gradient density 
field, (b) Frequency spectrum Sz{co,t = SOTq). The quantity Sq is defined as the time 
average of the main component So = (Sz(<x'0/ 0)t- 



order to allow a better observation of these three frequencies, we present on Fig. |3|b) a 
vertical cut of the time-frequency spectrum at time f/Tg = 50. As can be noticed, the 
frequencies satisfy the condition toi + CV2 = i^o- We will come back to this feature in 
the following sections. It is interesting to notice that although the two secondary fre- 
quencies seem to shift slightly in the long time range, the resonance condition prevails. 
In addition, the amplitude of the two secondary waves seems to decrease after some 
time. Two other frequencies, cv/N = and co/N = 0.98, have non vanishing contribu- 
tions to the signal. The first one corresponds to the mean flow generated by the plates 
of the wave generator via an Archimedes' screw type of entrainment. A second possible 
source of mean flow could be the reflection of the primary wave at the bottom of the 
tank. The other additional frequency, cv/N = 0.98, can be attributed to the non-lrnear 
interaction between the waves with frequencies L02/N = 0.24 and a;o/N=0.74. 

To extract more information on the various waves involved in this flow, one can filter 
the density gradient field around the three measured frequencies coq, loi and t02- This 
filtering operation is performed using the Hilbert transform method that was devel- 
oped in Mercier et al.| (2008|. This method consists in a first demodulation step with a 
Fourier transform in time, a time filtering around the desired frequency and then re- 
turning to real space by an inverse Fourier transform. In a second step, the signal is 
filtered in space, allowing separation of the waves corresponding to the four possible 
propagation directions for a given frequency. The result of this operation applied to our 
measured density fields is shown in Fig.|4] The filtered density gradient fields associated 
with each frequency are shown in the top row, at a fixed time. Note that for the three 
different cases, we have kept only one quarter of the possible wavevectors as shown 
in the sketches of Fig. |4] Phase and group velocities being orthogonal, the propagation 
direction is deduced by a rotation of 90°, the sense of rotation being chosen to have the 
vertical components of the phase and group velocities of opposite signs. For the left and 
center columns, it corresponds to a propagation from the the top left corner to the bot- 
tom right one, while for the right column the wave beam goes from the bottom right 
corner to the top left one. One can clearly observe that the filtering operation extracted 
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Figure 4: (color online) Top row: primary (left) and two secondary waves (center and 
right) obtained by applying the Hilbert spatial and temporal filtering to the density 
gradient field at f/To = 50. The sketches on top indicate the domain of wavevector 
{£,m) which are kept during the spatial Hilbert filtering for each column. Bottom row, 
corresponding phases associated to each frequency. The phase is displayed only where 
the wave amplitude is larger than 3% of the maximum. The experimental parameters 
are N = 0.91, cvq = 0.74 N and a = 0.5 cm. 



three distinct waves, each associated with its corresponding propagation angle. As men- 
tioned, the Hilbert transform allows identification of the orientation of the propagation 
of each wave, showing that the wave at frequency t02 propagates from right to left, con- 
trary to the two others. This explains why we observe signal outside the primary wave 
beam in the top left corner of Fig.|2|b), as mentioned in the previous section. 

This procedure also allows extraction of the phase ^/ of the signal at a given frequency, 
(pi{t,x,z) = cOjt zL ijX zLmjZ, where are the horizontal and vertical components of 

the wavevector fc, . This is shown at f = 50 Tq in the bottom row of Fig.ji] Again, one can 
clearly observe a pattern of stripes parallel to the direction of propagation of each wave, 
corresponding to a phase propagating in the perpendicular direction. At a fixed time 
and X (respectively z), the phase is linear with the position z (resp. x). The components 
ij and m, of the wave vectors for each wave can then be obtained by differentiating 
(pj{t,x,z) with respect to z (resp. x). For the experiment presented in Fig.|2j one obtains 
{h+h)/^o = 0.89 ± 0.17 and {mi+m2)/mo = 0.99 ± 0.07. It is striking that the three 

vectors also satisfy a resonance condition: kg = ki + k2 within the experimental uncer- 
tainties. Note that the larger error on the measurement of £, compared to m, is related 
to the fact that the secondary waves are rather horizontal. Consequently the horizontal 
projection of the wave beam, for LO2, gives us only one or two horizontal wavelengths. 
It yields a poorer measurement of the horizontal component of their wavevector. 

As a conclusion of this first set of experimental observations, we can say that the two 
secondary waves that are generated from the primary wave are the result of a resonant 
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triad interaction. We will therefore study analytically in the next section the conditions 
under which such a triad interaction can develop. 



3. Parametric subharmonic instability theory 

3.1. Derivation of the equations 

The two-dimensional d5mamics (in the x, z coordinates) of a Boussinesq fluid is usually 
described by the following two equations 

%*+H^H.,f) = '^^+.A'f, (3.2) 

where ip is the streamf unction, while b = gp/ p is the buoyancy perturbation and / the 
jacobian term defined as /(/i,/2) = ^xfi^zfi — ^zfi^xfi- The velocity field is expressed 
as "zT' = {-dztp,0,dx-ip). 

We seek solutions of the form 

b=f^Rjit)e'(^r~''i*'^+c.c., (3.3) 
j=0 

2 

= E 'i'j{t)e'^''r-'^'i^) + ex. . (3.4) 
/=0 

Introducing these solutions into Eq. | |3.1| and Eq. p.2\ , one gets 

2 

[Rj - iLOjRj + iN^ejYj]e'^''r^-''i^^ + c.c. = -}{b, xp) , (3.5) 

2 

\-k]{^^ - icvjl'j) - iljRj - VK^^ ■y^^y'-'^i'^ + ex. = -JiV^ip, ip) , (3.6) 

;=0 

if R denotes the derivative of the amplitude R. 

The usual inviscid linear dynamics of | |3.5| provides the polarization expression 

Rj = '-Yj for 7 = 0, 1 or 2 (3.7) 

with the dispersion relation 

Oj = SjN\£j\/^£j + mj (3.8) 



CO: 



where = ±1 defines the sign of the wave /. However, the linear system is forced on 
the right-hand side of Eqs. 13.5 1 and 13.6 1. After some manipulations, the Jacobian terms 
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can be written as 



+ c.c. 



(3.9) 



-[{ipniq - mp£,)K2YpY*]e'[(^p-^''>^-K-'^'))fl + c.c. . (3.10) 

We obtain now the evolution of a particular wavenumber component {k,-, co,-) associated 
with the stream function Y,-, in which r = 0, 1 or 2, by averaging both the left hand side 
and the right hand side over the period of that wave. The resonant terms on the right 
hand side that will balance the left hand side will be the waves fulfilling two conditions: 
a spatial resonance condition 



and a temporal resonance one 



h+k2 



(3.11) 
(3.12) 



Collecting resonant terms and using the polarization expression 1 3.7 1, the Jacobian 
term 1 3.9 1 can then be written as 



}{b,ip) = -(£iOT2-mi^2)N^ 

+ (£o'«2-OTo£2)N2 

+NRT, 



) YiY2e''*^'''''~'^''''' 



(3.13) 



where NRT stands for non resonant terms that are not important in the problem. In the 
same way, one gets the Jacobian term 1 3.10| 

/( V^ip, ip) = (£iOT2 - wxi^2) (k? - k|) Yi Y2e'(^o-^-'^o') 
-(£o'«2 - mo4)(Kg - K2)YoY^e'(^i-?— If) 

+NRT. (3.14) 

Introducing this result into equation | |3.6) , one obtains the three following relations be- 
tween Y,. and Rr for each individual phase exp[f (A:,- ■ r — 0Jrt)\ in which r = 0, 1 or 2 



(3.15) 
(3.16) 
(3.17) 
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where 70 = +1, 71,2 = — 1 and = [ipm^ — m^lq) (k^ — k^) , in which p,q,r = 0, 1, 2 or 
any circular permutation. 

3.2. Slow amplitude variation 

Experimental results have strongly suggested that the amplitude of T varies slowly 
with respect to the period of the primary wave. It is therefore appropriate to consider 
< fwy Yy so that differentiating (jsi/j leads to 

07; 



(3.18) 



One gets for Eq. ||3.5 



{N^el - cvIkDWo - 2icvo4% - icvQVK^Yo + /a;oSoao(YiT2) 



{£im2 - mi£2)N^ 



{N^£^ - CVrKr)^r - 2za;,<Y, - icVrVK'^Yr + fa;,Sra,(YoY*) 

Each wave satisfies the dispersion relation, N^£-^ = cv}k}. After some calculations, one 
gets 



-{£omp-moep)N^ 



(3.19) 



(3.20) 



2za;oK^Yo 

= i{£im2 — m\(.2) 
which can be simplified as 



icOQSQOLQ r(^l7«2 — mii2)N 

I \ COl ^2 

soa;o(K?-K|) + 4N2(^-^ 

COl ^2 



T1T2 - fa;oVK*Yo (3.21) 
YiY2-fa;oVK^To,(3.22) 



To = IoYiY2 - -K^To. 



Similar calculations for waves 1 and 2 lead to 



Y2 = -l2YoTj--K^Y2, 



where 



Ir = 



IcOrK} 



COr{Kl-K^^) + lrN^ 



3.3. Solution 



(3.23) 

(3.24) 

(3.25) 

(3.26) 



We consider that Yq corresponds to the primary wave and is constant in early times 
since amplitudes of the secondary waves, Yi and Y2 are negligible with respect to Yq. 
One can combine equations ( 3.24^ and 1 3.25| to get 

2 

Yi = hh^l^i - ^K^KfYi -^{4 + k|)Yi . (3.27) 



The solution of Eq. 13.27 1 leads to the expression Yi^2(T) = ^1,2 + Bi,2e''- , by 
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Figure 5: (color online) The curves represent the location of {^i, nii) satisfying Eq. 13.29 

for the three possible combinations of signs, once the wavevector = {(-a, mo) of the 
primary wave is given. The curves are represented by solid lines when the growth rates 
have a positive real part, while the dashed lines corresponds to the neutrally stable 
cases, for which the real part of the growth rate is always zero in a non viscous case and 

negative when viscosity sets in. Two examples of vector triads (fcg , ki,k2) are shown. 



introducing the growth rate 

^± = -l{4 + 4) ^\1y6^A- 4Y + hh\^o?. (3.28) 

A vanishingly small amplitude noise induces the growth of two secondary waves. In 
conclusion, a primary plane wave can be unstable by a parametric subharmonic mech- 
anism. The growth rate of the instability is not only function of the characteristics of the 
primary wave, namely its wavevector and its frequency through the coefficients li and 
I2, but also its amplitude Tq and the viscosity v. 

3.4. Resonance loci and growth rates 

In the following, we consider that only the primary wave Yq, of given frequency coq, 

wavevector fcg = (^O/'^o) ^rid sign sq, is present initially in the system, while Y12 is 
experimentally only present as the noise. Without loss of generality, we can choose 
So = +1 (which is given by the experimental configuration, or more precisely by the 

direction of rotation of the cams). The two secondary waves (si, coi, ki ) and (s2, 0J2, k2) 
which could form a resonant triad with the primary wave have to be determined us- 
ing the resonance conditions | 3.11| and | 3.12| . From the dispersion relation for internal 



waves 13.8 1, the resonance conditions lead to 



So — , = Si — + S2 , — ^ = ■ (3.29) 

For a given primary wave (sq, io, mo), the solution of this equation for each sign com- 
bination (so,si,S2) is a curve in the {ii,mi) plane which is presented in Fig. p] Since 
Sq — +1, it is then necessary to consider four sign combinations for (si,S2): (— , — ), 
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Figure 6: (color online) Growth rates a computed from Eq. 1 3.28 1, as a function of the 
wave vector modulus kj in panel (a) and as a function of the wave frequency coi in panel 
(b). The three possible sign combinations have been superimposed. The primary wave 
vector kf) is chosen arbitrarily, as well as its amplitude, which value, normalized by the 
viscosity, is Yq/v = 100. Kg and loq are respectively the wave vector modulus and the 
frequency of the primary wave. 



(+,—), (— , +) and (+, +). However, there is no solution of equation 13.29 1 in the (— , — ) 
case. In addition, the combinations (+, — ) and (— , +) are neutrally stable, i.e. the real 
part of cr is always zero in a non viscous case and negative when viscosity sets in. For 
this reason, in what follows, we will focus only on the (+, +) combination. 

In Fig. [sj any point of the solid curve corresponds to the tip of the ki vector, wich 



satisfy the equation [3.29 with sg = = S2 = +1. We obtain also k2 by construction, 
closing the triangle. Two examples are illustrated with solid arrows. However, to predict 
which is the one expected to be seen experimentally, one has to determine the largest 
growth rate. For a given sign of m-[, three distinct parts of the (+, +) curves can be 
observed in Fig. |5| corresponding to ii/£o > 1, £i/£o < and < ii/io < 1- It will 
prove convenient in what follows to separate the study in these three regions, in order 
to evidence the direction of energy transfer (towards smaller scales or larger scales). 

In Fig. |6| the predicted growth rates a are plotted as a function of / kq and coi / o^q 
only for the (+,+) combination, with distinct line types for the different regions defined 
above. There are two curves for each line type, corresponding to different signs of mi 
(top half and bottom half in Fig.|5|. 

One can notice that in both plots, the solid curve {£i/£o > 1) ai^d the dashed curve 
{£l / £o < 0) reach the same maximum value. Moreover, one can observe in Fig.|6|b) that 
the maximum are obtained for coi = 0.37 loq (dashed curve) and for coi — 0.63 ojq (solid 
curve), two frequencies which sum is equal to Wq. More generally, Fig.|6p3) shows that 
the solid and dashed curves are mutually symmetric with respect to the central value 

wl wq = 0.5. It shows that if k\ is selected on the solid curve, ki will be selected on the 
dashed curve and conversely. In an arbitrary way, we will always take k\ on the solid 

curve (largest wavenumber) and consequently ki will be determined by the dashed 
curve (lowest wavenumber among the secondary waves). 

As to the dotted dashed curve, one can notice that it has two maxima. The two sec- 
ondary wave vectors are, in this case, selected in the same area (0 < l\l < 1). For all 
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the curves, one notices that the growth rates become negative when ki ^ oo (because 
of viscosity). 

Interestingly, the growth rate is positive for a broad range of wavenumbers and it 
is rather flat on the maximum, emphasizing that the parametric resonance is weakly 
selective in this regime. The values of Ki corresponding to significant growth rates are 
of the same order of magnitude as the primary wavenumber kq, indicating that the 
viscosity has a significant effect on the selection of the excited resonant triad, preventing 
any large wave number secondary wave to grow from the instability. For the frequency 
value considered, the maximum growth rate is obtained for ki =2.5 kq and K2 = 1-62 kq. 



4. Comparison with the experiment 

4.1. Wave vectors 

In section |2j we observed that the experimental values of frequencies and wave vectors 
satisfy the conditions of temporal and spatial resonance | |3.11[ and 1 3.12 1. The theoretical 
derivation gives us the dependence of the growth rate with the wave vector of one of 
the secondary waves (the second wave is then defined by the resonance condition). We 
can now check if the observed wave vectors correspond to the largest growth rates, as 
could be expected. For pedagogical reasons, the value of Yq used to compute the plots in 
section|3]was arbitrarily chosen to enhance the case (external branches of the resonance 

loci) where the energy transfer goes from large scales to smaller scales. The vector Rq 
was also chosen arbitrarily. In this section, we will use the experimental values for fcg 
as well as for the amplitude Yq. This last quantity is measured using the modulus of 
the Hilbert Transform of the density gradient field for the primary wave, calculated in 
the same area used for the measurement of the growth rate. Figure [Zja) shows the the- 
oretical resonance loci (blue curve) computed identically to the one in Fig.|5jbut usme 
the experimental kg . The three arrows correspond to the experimental measurement of 
the three wave vectors. As mentioned in section |2j it can be observed that the spatial 
resonance condition is satisfied. Figure [7p3) represents the evolution of the growth rate 
with Ki. Comparing Fig.|6|a) and [Tjb), is interesting to note that the reduction of the pri- 
mary wave amplitude by a factor of 4 results in an enhancement of the (0 < l\l Iq < 1) 
type of instability, with respect to the {ii/io > 1) type. 



According to the theory, the secondary waves generated by the instability should be 
the ones corresponding to the largest growth rate. However, an instability is initially 
generated via a very small amplitude noise and because of the specific experimental 
conditions, not all temporal and spatial frequencies will be present in this noise with 
the same amplitude. In addition, as can be seen in Fig. [Tjb), the growth rate curve in 
the case of the (0 < £i/£o < 1) type of instability, although it has a slightly higher max- 
imum, is much less broad than the curve for the other type, therefore showing a very 
strong wavenumber sensitivity. For these reasons, we will introduce a less constrained 
selection criteria by allowing the unstable wave vectors to be associated with a growth 
rate belonging to an interval between 90% and 100% of the maximum growth rate. This 
selection criteria is illustrated in Fig.|7](b) by thicker line and the loci for the correspond- 
ing wave vectors are represented in Fig.[7|(a) by thick black segments. It is interesting 
to note that this procedure allows a selection of wave vectors in different regions of the 
resonance loci curve. The measured wave vectors, within their experimental error, fall 
on one pair of these selected regions, namely on the external branches of the resonance 




Figure 7: (color online^ (a) The solid blue curve represents the resonance loci for the sec- 
ondary wave vector ki , in the case where the primary wave vector fcg corresponds to the 
experiment presented in §|2| i.e. = 66 ± 2 m^^, mq = 59 ±2 m^^ and cvq/N = 0.74. 
The three arrows represent the experimental measurement of the three wa ve ve ctors 
k() (black), ki (red) and k2 (green), (b) Growth rates a computed from Eq. 1 3.28 1, as a 

function of the wave vector modulus, in the case where the primary wave vector kg 
corresponds to the experiment presented in section |2j as well as the primary wave am- 
plitude (Tg/v = 28). The two horizontal lines mark the maximum growth rate and 90% 
of the maximum. The thick lines on the left panel show the regions corresponding to 
a growth rate belonging to this interval. The experimentally measured growth rate are 
shown for the observed wavenumber ki (red square) and K2 (green circle). The experi- 
mental error bars on the growth rate measurement beerng large (of the order of 25%, see 
text), they are not represented for these two points, for the sake of clarity. 



loci. This selection corresponds, as mentioned earlier, to a slightly lower growth rate 
than the maximum, but to a broader curve. 

4.2. Growth rates 

We will now compare the experimental and theoretical values of the growth rate. The 
experimental determination of the growth rate a is obtained by calculating a time- 
frequency spectrum of an area taken inside the primary wave beam, as defined in Fig. [l] 
by the grey area. Then, the time evolution of each partial component frequency is plot- 
ted. Figure|8]shows this evolution for the same experimental run used in § 2] The growth 
rate is then defined as the slope of the growth region. It must be noted that the ampli- 
tude of the primary plane wave is not constant on the entire beam because of the viscous 
damping. For this reason, the position of the area where the spectrum is computed will 
influence the result for the determination of a. An estimate of the error due to this effect 
is therefore performed by varying this position. 



For the theoretical evaluation of the growth rate, a necessary ingredient is the ampli- 
tude Yq of the primary wave. This amplitude cannot be obtained directly from the am- 
plitude of oscillation of the generator plates, since there is a conversion factor between 
both described by |Mercier et al.| ( 2010| . This factor is not well known, and strongly de- 
pends on the wave frequency, because of the change of angle it implies of the wave 
propagation, while the generator motion is only vertical. For this reason, Yg has to be 
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Figure 8: (color online) Amplitude of time frequency representation, Sz{oJ,t) of the sec- 
ondary wave coi and 0J2 normalized by the amplitude of the primary wave. The dashed- 
dotted lines show the linear fit, from which the experimental value of the growth rate 
a is extracted. The quantity So is defined as the time average of the main component 

So = {Sz{cvo,t))t- 

measured experimentally, introducing an experimental error in the prediction. It leads 
to a theoretical prediction for the growth rate cr = 0.087 ± 0.005 s^^, while the exper- 
imental measurement gives a = 0.08 ± 0.02 s~^. The agreement between theory and 
experiment is rather good. 

5. Dependence with primary wave frequency and amplitude 

So far the experimental results shown correspond to a given amplitude, frequency 
and wave vector for the primary wave. It seems appropriate to study the effect of a 
variation of these parameters. 

First, we have observed, based on the theoretical study, that the value of kq has no 
effect on the selection of the secondary waves. It leaves unchanged the shape of the 
growth rate evolution with kj and cvi (Fig. |6}. The only effect of an increase of kq is an 
increase of the overall value of the growth rate. For this reason, it was not necessary to 
perform a systematic experimental study of the dependence of the growth rate with kq. 

In a first part, we will discuss the dependence with the primary wave frequency a;o- 
Figure |9] shows the evolution with cvq/N of the (normalized) norms Ki and K2 of the 
secondary waves corresponding to the maximum growth rate. One pair of curves (solid 
lines) corresponds to wave vectors on the external branches of the resonance loci and 
the other pair (dashed-dotted lines) to wave vectors on the central part of the resonance 
loci. These curves are obtained in the case where the primary wave amplitude is fixed at 
the value of the experimental run described in section|2] We can observe that in this case, 
when the secondary wave vectors belong to the external branches, there exist a range of 
frequency (for coq/N > 0.78) where the secondary wavelengths are both smaller than 
the primary wavelength, corresponding to an energy transfer to smaller scales. On the 
contrary, when the secondary wave vectors belong to the central part of the resonance 
loci (dashed-dotted lines), one of the secondary wavelengths is always larger than the 
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Figure 9: Evolution of kj / kq and K2 1 Kg with the pulsation a;o for a fixed amplitude. Solid 
curve: Configuration where the triad is formed on the external branches of the resonance 
loci (see illustration in the top inset). Dashed dotted curve: configuration where the triad 
is formed on the central part of the resonance loci (see illustration in the bottom inset). 
The thick part of curves represent the region where the growth rate is larger in this 
configuration, compared to the other. 



primary wavelength, while the other is smaller, allowing for energy transfer to both 
smaller and larger scales. 

Out of these two sets of secondary wave vectors, the one which corresponds to the 
largest growth rate is illustrated on the graph by a thicker line. There is therefore a 
transition frequency (in this case around cvq/N = 0.77) for which the instability "jumps" 
from one solution to the other. 

As mentioned, this change in behavior is observed for a given value of the ampli- 
tude Yq. As this amplitude varies, the transition frequency changes. Figure 10 shows 
(solid line) the evolution of this transition frequency as a function of the amplitude. 
Above this curve, the secondary waves are selected on the outer region of the resonance 
loci, while below the curve the secondary waves are selected in the central zone. 

On the same figure, we superimposed the experimental points obtained for different 
values of loq and Yq, showing in which cases the parametric subharmonic instability 
(PSI) was observed. Three different wave generator amplitudes were used to produce 
these series of points (cam eccentricities of 0.1, 0.5 and 0.75 cm). These three series corre- 
spond to the three almost vertical series of black symbols on the graph. One must keep 
in mind that the wave amplitude is computed from the measured signals, therefore a 
fixed amplitude of the wave generator does not correspond exactly to a fixed wave am- 
plitude. 



It can be noted that for large enough excitation amplitudes, PSI is observed, but 
mostly above the transition line. Actually, after careful observation, it turns out that 
the two points with PSI below this line are also case of PSI with wave vectors on the 
external branches of the resonance loci. This feature may be related to the uncertainty 
on the frequency selection mentioned in section 4.1 
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Figure 10: (color online) Phase diagram showing the dependence of the PSI instability 
with excitation frequency and primary wave amplitude. The solid line represents the 
transition between the configuration where the triad forms on the external branches 
of the resonance loci and the configuration where the triad is on the central part. The 
dashed line represents the cutoff frequency under which energy transfer to larger scales 
is involved in the PSI. Symbols correspond to experimental data points: (o) observed 
PSI, (+) no apparent PSI. 



Our guess is that PSI cannot develop in the present experimental tank when one of 
the secondary waves has a larger wavelength than the primary wave. No cascade to 
larger scale is allowed. The reason for this behavior could be the too small number of 
wavelengths which is excited by the generator. 

For the smallest excitation amplitude, no PSI was observed during the recording, even 
for experimental runs located above the solid line of figure 10 Equation | 3.28| shows 
that unstable solutions are possible only if the amplitude of the primary wave exceeds a 
viscous threshold value given by 2 ^vk\Ki/ (hh)- However, the excitation amplitude is 
one or two orders of magnitude larger than this viscous threshold which cannot explain 
the absence of the instability. This behavior could again be related to the size of the ex- 
pected wavelengths at this amplitude. Indeed, it can be observed in Fig.|9]that the lower 
solid line crosses the ordinate 1 at a given frequency. This implies that, below this fre- 
quency, although the secondary waves are selected on the external branches, a transfer 
to larger scales is however involved. We have computed theoretically, as a function of 
the wave amplitude, the threshold frequency under which a transfer to larger scales is 
present for secondary waves selected on the external branches. It is plotted in Fig. 10 as 



a dashed line. The small amplitude series, where no PSI is observed, happens to be on 
the left of this line, therefore in the region involving transfer to larger scales. As in the 
case of small frequencies, this might explain the absence of observed PSI. To summa- 
rize, the only region in the phase space where no energy is transfered to larger scales is 
the top right comer (above the solid line and right of the dashed line) and it is the only 
region where we observe PSI, which seems to agree with our assumption of limitation 
due to the finite size of the primary wave beam. 
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6. Discussion 

Using a synthetic schlieren technique in a linearly stratified fluid, we have observed 
the production, from a primary internal plane wave, of two secondary waves of smaller 
frequency and wavelength. The mechanism at play is a 3-wave resonant interaction, 
satisfying a temporal resonance condition for the three pulsations coq = co-i + 0J2 and a 

spatial resonance condition for the three wave vectors, ko = ki + k2- This resonant in- 
teraction can be developed analytically, in order to produce theoretical predictions that 
we compare to our measurements, which are derived from applying a time-frequency 
spectrum and a Hilbert transform to our data. A good agreement was found for the 
wave vectors of the secondary waves, as well as for the instability growth rate. 

The relevance of PSI for oceanic internal waves has been recently renewed by several 



numerical simulations with realistic conditions (see for example Hibiya et al. ( 20021 



or 



MacKinnon & Winters| ( |2005[ ), which have been followed by very mterestmg heli 



9 



measurements which also pointed to the presence of PSI activity in the ocean. One might 
in particular refer to measurements close to Hawaii by Rainville & Pinkel ( |2006a| . Carter 
' & Gregg| ( |2006| or |Alford etal.| ( |2007| |. 



However, in the ocean, due to the large spatial scales, the role of viscosity is reduced 
compared to the experiment. In order to illustrate this point, let us consider a typi- 
cal oceanic internal wave, measured in situ by observing the periodic oscillations of 
isotherms (see figure 1 in Cair ns & Williams | ( [T$7 6^I or figure 6.13 in |Sutherlan"d| ( |20T0) ). 



Theses oscillations have an amplitude of about 20 m, with a period of about half a day. 
Ocean internal waves t5^ically have wavelengths from hundreds of meters to tens of 
kilometers (Garrett & Munk 1972{ |LeBlond & Mysak|1978||Susanto etal.|2005[ |: as a typ- 



ical value, we shall take 1000 rnTCompared to our experimient, it yields an amplitude of 
the stream function Yq / v about four orders of magnitude larger. We therefore recom- 
puted as an example (see Fig. [TT| the curves of Fig. [Tjb), as well as the corresponding 
plots as a function of the wave frequency, for an amplitude 10000 times larger than the 



experiment. As can be observed in Fig. 11 a), the most unstable mode corresponds to 



the situation where the norm of both secondary wave vectors are much larger (10 to 20 

times in our example) than the norm of the primary wave vector, leading to k^ ~ k2 
and correspondingly (Fig. [TTpj)) to coi LO2 — coo/2. This is what led historically to 
call this resonant interaction "Parametric subharmonic instability", because of the sim- 
ilarity with parametric instability, like that occurring for a pendulum with oscillating 
suspension. 

In our experimental tank, with dimensions orders of magnitude smaller than in the 
ocean, we observe that viscosity has a much larger effect on the selection of the unstable 
modes. Indeed, because of viscosity, Ki and K2 are of the same order of magnitude as 
kq. In addition, a second effect of viscosity is to allow two different behaviors for the in- 
stability. First, the instability can generate two new internal waves whose wavelengths 
are smaller than the primary wavelength. We have then an energetic transfer to smaller 
scales, as is the case in the ocean. But theoretically, when Yq / 1/ is small enough, the in- 
stability can also generate two secondary internal waves with one of the wavelengths 
larger than the primary wavelength, and the other one smaller. In this case, the ener- 
getic transfer will be to larger and smaller scales simultaneously. However, we never 
observed experimentally this second behavior, although it was expected from the ana- 
lytical study. We observe the instability, on a plane wave, only when the theory predicts 
two smaller wavelengths. It is possible that the finite size of the beam impedes the de- 
velopment of the secondary wave, which has a larger wavelength. It could be interesting 
to perform some numerical simulations, with a confined beam to confirm this behavior. 
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Figure 11: (color online) Growth rates tr as a ftinction of the wave vector modulus 
in panel (a) and as a function of the wave frequency coi in panel (b) for typical oceanic 
parameters (see text). 

Another issue that needs to be discussed is the physical location of the birth of the 
secondary waves. In all the experiments performed, the instability appears first in a 
region very close to the wave generator. This is understandable because it is the location 
where the primary wave first appears and where its amplitude is maximum. Then it 
occupies the whole volume. 

In addition to the importance of the mechanism for the dissipation of waves, an effect 
on the background media in which waves are propagating can be expected: indeed, the 
primary wave could give rise through this instability to two secondary waves which 
are above the overturn threshold. Because of the latter, the stratification will evolve 
through these mixing events, influencing through this feedback effect the threshold for 
the mechanism itself. 
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